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ABSTRACT 

The dark energy that appears to produce the accelerating expansion of the universe 
can be characterized by an equation of state p = wp with w < —1/3. A number of 
observational tests have been proposed to study the value or redshift dependence of 
w, including SN la distances, the Sunyaev-Zel’dovich effect, cluster abundances, strong 
and weak gravitational lensing, galaxy and quasar clustering, galaxy ages, the Lya for¬ 
est, and cosmic microwave background anisotropies. The proposed observational tests 
based on these phenomena measure either the distance-redshift relation d(z), the Hub¬ 
ble parameter H(z), the age of the universe t(z), the linear growth factor D\(z), or 
some combination of these quantities. We compute the evolution of these four observ¬ 
ables, and of the combination H(z)d(z) that enters the Alcock-Paczyznski anisotropy 
test, in models with constant w, in quintessence models with some simple forms of the 
potential and in toy models that allow more radical time variations of w. Mea¬ 

surement of any of these quantities to precision of a few percent is generally sufficient 
to discriminate between w = — 1 and w = —2/3. However, the time-dependence pre¬ 
dicted in quintessence models is extremely difficult to discern because the quintessence 
component is dynamically unimportant at the redshifts where w departs substantially 
from its low -2 value. Even for the toy models that allow substantial changes in w at 
low redshift, there is always a constant-u; model that produces very similar evolution 
of all of the observables simultaneously. We conclude that measurement of the effective 
equation of state of the dark energy may be achieved by several independent routes in 
the next few years, but that detecting time-variation in this equation of state will prove 
very difficult except in specialized cases. 


Subject headings: cosmology: theory 


1. Introduction 

The big cosmological surprise of recent years is that the dominant form of energy in the 
universe has negative pressure and is therefore causing the expansion of the universe to accelerate. 
The most direct evidence for acceleration comes from the Hubble diagram of Type la supernovae 
(SN la), in particular the relative apparent brightness of SN la at redshifts z ~ 0 and z ~ 0.5 — 1 
(Riess et al. 1998; Perlmutter et al. 1999). However, other strong arguments for a “dark energy” 
component follow from combining the cosmic microwave background (CMB) evidence for a spatially 
flat universe (Netterfield et al. 2001; Pryke et al. 2001) with either a minimum age to ~ 13 Gyr 
(Vandenberg, Stetson, & Bolte 1996) or dynamical evidence that the density of clustered matter 
is well below the critical density (see Bahcall, Fan, & Cen 1997; Carlberg, Yee, & Ellingson 1997; 
Weinberg et al. 1999b for examples of three distinct routes to this conclusion, though there are 
many others). The first combination, together with a Hubble constant Ho ~ 70 kms -1 Mpc -1 = 
(14 Gyr) -1 (Freedman et al. 2001), requires a component whose gravitational acceleration roughly 
cancels the gravitational deceleration caused by the pressureless matter, so that to ~ Hjj 1 . The 
second combination requires that the dominant form of energy be unclustered, though it implies 
nothing more specific about its equation of state. A more model-dependent argument for a negative 
pressure component comes from the success of inflationary models with cold dark matter (CDM) 
and a cosmological constant (A) in matching a variety of constraints from CMB anisotropies and 
large scale structure measurements (see Wang, Tegmark, &; Zaldarriaga 2001 for a recent review). 

In this paper, we explore the prospects for determining the equation of state of the dark 
energy component through a variety of observational methods. A true cosmological constant can 
be treated as a vacuum energy with time-independent density and pressure related by p = —p. 
Current observations favor an equation of state fairly close to this prediction (Garnavich et al. 
1998). However, a number of authors have considered the more general possibility that the negative 
pressure component is a scalar field (a.k.a. “quintessence”) with energy density determined by 
its potential and effective equation of state p = wp, where w can be constant or time-varying 
(Ratra & Peebles 1988; Turner & White 1997; Caldwell, Dave, & Steinhardt 1998). Interest in 
models with time-varying w has been spurred by arguments that certain simple potentials lead 
“naturally” to a negative pressure quintessence component that dominates the expansion at late 
times, independent of the initial conditions (Zlatev, Wang, & Steinhardt 1999; Steinhardt, Wang, Sz 
Zlatev 1999). Variants on this theme include fields with a non-standard kinetic term (Armendariz- 
Picon, Mukhanov, & Steinhardt 2000) or models with a complex scalar field (Boyle, Caldwell, & 
Kamionkowski 2001). 

Further afield, there is the possibility that the negative pressure component is a network of 
frustrated topological defects (Vilenkin 1985; Spergel & Pen 1997), or that cosmic acceleration 
arises from a breakdown of general relativity rather than the addition of a new energy component 
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(Mannheim 2001; see also Tegmark 2001). The hope, thus far unrealized, is that one of these ideas 
will eventually provide a natural explanation of why the vacuum energy density is 120 orders-of- 
magnitude below the Planck scale and why it is comparable to the matter density at the present 
day, without having to resort to anthropic selection arguments (Efstathiou 1995; Martel, Shapiro, 
& Weinberg 1998). 

Any clear evidence that w ^ — 1, or, better still, that w varies in time, would provide crucial 
clues towards understanding the physics of the dark energy. Through its influence on the cosmic ex¬ 
pansion history, this component affects many observable phenomena, including CMB anisotropies, 
the Lya forest, strong and weak gravitational lensing, the anisotropy of quasar and galaxy cluster¬ 
ing in redshift space, the ages of the oldest galaxies as a function of redshift, and standard-candle or 
standard-ruler measurements of the distance-redshift relation. This paper discusses these potential 
observational tests in a unified fashion. The equation of state determines the history of the energy 
density which, together with the densities p rn and p r of matter and radiation, determines the 
evolution of the Hubble parameter H(z) via the Friedmann equation. The history of H(z) in turn 
determines the age of the universe t(z), the growth factor of linear perturbations D\(z), and dis¬ 
tance measures like the angular diameter distance cIa{z) or luminosity distance cIl(z), which are 
related to each other by cosmology-independent powers of (1 + z). Essentially all proposed tests of 
the properties of the negative pressure component amount to measurements of H(z ), t(z), D\(z), 
or d(z), or some combination of them, at redshifts accessible to a particular observational tech¬ 
nique. We will investigate the dependence of these four quantities, and of the specific combination 
H(z)cIa{z ) that is constrained by the Alcock-Paczynski (1979; hereafter AP) anisotropy test, on 
the value and time history of w. 

Our paper joins, and, we hope, complements, a flood of recent papers that examine the 
prospects for specific tests and specific data sets in much greater detail. Since the strongest 
evidence for A or a quintessence component comes from SN la observations, and substantial im¬ 
provements are likely from ground-based campaigns and possibly a dedicated satellite (SNAP; see 
http://snap.lbl.gov), many authors have examined the extent to which present or future SN la 
observations can constrain w(z) (Turner & White 1997; Garnavich et al. 1998; Astier 2000; Chiba 
& Nakamura 2000; Huterer & Turner 2001; Saini et al. 2000; Barger & Marfatia 2001; Chevallier & 
Polarski 2001; Maor, Brustein, & Steinhardt 2001; Ng & Wiltshire 2001; Podariu, Nugent, & Ratra 
2001; Wang & Garnavich 2001; Wang & Lovelace 2001; Weller & Albrecht 2001). Because CMB 
anisotropy predictions depend most strongly on the sum of p$ and p rn while SN la distances depend 
more nearly on the difference, the combination of these complementary observations yields much 
tighter constraints on the negative pressure component than either does alone (Caldwell, Dave, 
& Steinhardt 1998; Efstathiou 1999; Baccigalupi et al. 2001; Corasiniti Sz Copeland 2001; Doran, 
Lilley, & Wetterich 2001). The Sunyaev-Zel’dovich effect or size of radio sources offer alternative 
ways of measuring (Ia{z) (Birkinshaw 1999; Lima & Alcaniz 2001), and the volume-redshift test 
using galaxy counts constrains the combination d\(z)H~ l (z) (Newman & Davis 2000, 2001). The 
evolution of the galaxy cluster mass function can constrain the linear growth factor Di(z) (Benabed 
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& Bernardeau 2001; Doran, Schwindt, & Wetterich 2001; Haiman, Mohr, & Holder 2001; Newman 
et al. 2001; Weller, Battye, & Kneissl 2001), and population synthesis modeling of galaxy spectra 
can constrain t(z ) (Lima & Alcaniz 2000). Jimenez &; Loeb (2001) suggest that relative galaxy 
ages can be used to measure dz/dt , and thus H(z). Hui (1999) and Huterer (2001) have examined 
constraints on w that can be obtained from weak lensing, while Hu (2002) has considered lensing in 
combination with the CMB. Calvao, de Mello Neto, & Waga (2001) have discussed constraints that 
could be obtained by applying the AP test to the 2dF quasar redshift survey of Boyle et al. (2000; 
for related discussions see Hui, Stebbins, & Buries 1999; Cappi 2001; Dalai et al. 2001; McDonald 
2001 ). 

Most of these papers have considered the potential observational constraints singly, or in pairs. 
The goals of our more abstract discussion, where we consider all of these observables together but 
do not focus on specific observational strategies, are twofold. First, we aim to understand what level 
of precision is necessary with any of these quantities to obtain useful constraints on w. Second, 
we want to know whether these different observables provide complementary information about 
the time-variation of w, breaking degeneracies that exist for a single measure by probing different 
aspects of the expansion history. Unfortunately, our conclusions on the latter point are pessimistic 
— there are many different ways to measure w, but distinguishing a time-varying w from a constant 
w is likely to prove difficult. The papers by Wang et al. (2000) and Tegmark (2001) also consider 
multiple observables, focusing on present constraints and future prospects, respectively. Tegmark’s 
paper, in particular, is similar in spirit to ours, but different in the way that it frames the problem 
and evaluates the prospects. 

In the next section we discuss the various quintessence models that we examine in this paper. 
We discuss the observables in §3, beginning with the formulas that relate these quantities to the 
expansion history and proceeding to a brief account of observations that might measure these 
quantities in the next few years. We present our results in §4, first for the quintessence models 
described in §2, then for a class of “toy” models designed to allow stronger time-variation of w at 
low redshift. We summarize our conclusions in §5. 


2. Quintessence Models 


We will adopt the language and calculational framework of quintessence models, though most 
of our general conclusions are also relevant to other possible explanations of cosmic acceleration, 
like those mentioned in §1. Also, in light of evidence from the location of the first acoustic peak 
in the CMB anisotropy spectrum (Netterfield et al. 2001; Pryke et al. 2001), we will restrict our 
attention to spatially flat models. 


The Friedmann equation for a spatially flat, expanding universe can be written 


- = H(z) = H 0 


Pi(z) 
Pi, 0 


( 1 ) 
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Here a is the scale factor, a is the derivative of the scale factor with respect to time t, Hq is the 
value of the Hubble parameter at the present time to, and is the present density of some zth 
component of the energy density relative to the present critical density (fi^o = Pi,o/p c ,o)- For 
adiabatic expansion, the energy density of a component with equation of state p t = wipi with 
constant Wi evolves with redshift as 

^ = (1 + z) ni , (2) 

Pi, 0 

rii = 3(1 + Wi). (3) 

Normal matter has Wi = 0 and n t = 3, while radiation has Wi = 1/3 and n* = 4. A true cosmological 
constant, with pi =const, n t = 0, has Wi = —1. We will often refer to models in terms of the energy 
density scaling index n, defined by equation (2), rather than by w itself, since the value of n more 
directly captures the impact of a component on the expansion history. 

A coasting expansion, in which comoving observers have constant velocity, has H(z) oc (1 + z). 
An accelerated expansion requires, at a minimum, that the dominant energy component have n* < 2, 
and thus Wi < —1/3. [More precisely, (w), the density-weighted average value of w, must satisfy 
(w) < —1/3]. Quintessence, a term reintroduced to cosmology by Caldwell, Dave, & Steinhardt 
(1998) after millennia of neglect, refers generically to a scalar field with equation of state p$ = 
and < 0. The first class of models that we consider are those in which w$ is constant. In this 
case, the Friedmann equation can be written 

H(z) = — = Hq o(l + z) 4 + fl m o(l + z ) 3 + fij, o(l + z) n , (4) 

a v 

with n given by equation (3). 

More general models often treat quintessence as a minimally coupled scalar field (/), obeying 
the equation 


where w for the scalar field is 


= -3 H<j> - 


( 1 / 2 )4> 2 - V{(t>) 


( 1/2 )<p + V(<!>)' 

When V (< j>) is an exponential or a negative power-law, the scalar field has the desirable property that 
its final evolution is independent of initial conditions, a behavior that has been dubbed “tracking” 
(Zlatev, Wang, & Steinhardt 1999; Steinhardt, Wang, &; Zlatev 1999). The negative power-law 
potentials lead to constant w when the contribution from the scalar held energy density is sub¬ 
dominant (Ratra & Peebles 1988; Liddle & Scherrer 1999), but when the scalar held energy density 
comes to dominate at late times, the value of w changes. In principle, then, such models should be 
observationally distinguishable from models with constant w. 

For our second class of models, we have chosen a subset of the negative power-law potentials, 
where 

R( 0 ) oc r, (7) 
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Fig. 1.— The evolution of n = d\og p^/d\og(l + z) as a function of redshift z, for a scalar field 
with the indicated potential and a cosmological model with = 0.4, Q r _o = 9.8 x 10 5 , and 

= 1 — 0 — fir,o ~ 0.6. In this figure, and in all other figures in the paper, the horizontal 

axis effectively represents a set of discrete bins, so it is neither linear nor logarithmic. 

with a < 0. If the dominant component has a density that scales as p oc (1 +z) m (e.g., m = 4 during 
the radiation-dominated era and m = 3 during the matter dominated era), then these models have 

n=[a/(a — 2)\m (8) 

when ptf, <C p m (Liddle & Scherrer 1999). At late times, when the scalar field energy density begins 
to dominate, equation (8) no longer holds, and n changes with time. We have chosen to examine 
two representative cases: a = — 1 and a = —6. This choice is somewhat arbitrary, but these are 
the same cases that are discussed by Zlatev, Wang, &; Steinhardt (1999). 

The evolution of n in these models is displayed in Figure 1, assuming cosmological parameter 
values n m>0 = 0.4, ll r .o = 9.8 x 10 -5 , and = 1 — flm.o ~ fl r ,o- (This value of Q r ,o corresponds 
to a photon temperature of T = 2.73 K, a standard neutrino population, and a Hubble parameter 
of H 0 = 65 kms -1 Mpc -1 . This is the only place in our calculations where Ho enters, and it has 
a very small effect on our results). For these models, we define n(z) to be the local logarithmic 
derivative of p<f, with respect to (1 + z). Figure 1 shows that the value of n for 3 < z < 10 is almost 
exactly constant and given by equation (8), namely n = 1 for a = —1 and n = 9/4 for a = —6. At 
z < 3, n decreases slightly, reaching present-day values of n = 0.77 for a = — 1 and n = 1.89 for 
a = —6. In §4 we will see whether cosmological tests can detect these slight changes in n. 


Although models with an exponential potential are quite natural and can lead to p$ ~ p m at 
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all times, they are ruled out for several reasons: they tend to give p^ too large during primordial 
nucleosynthesis (Ferreira & Joyce 1997), and they lead to n = 3 at late times, which does not 
produce an accelerated expansion. These problems are remedied in the model of Albrecht & Skordis 
(2000), who introduced a potential consisting of an exponential multiplied by a polynomial, 

V(ct>) = [(<f>-By+ A]e- X ‘t’, (9) 

where A, B, 7, and A are constants. These constants can be chosen to produce a model for which 
n = m at early times, when the scalar field is sliding down the exponential potential, but n = 0 
at late times, when the scalar field settles into the local minimum in the potential. The constants 
in this model must still be tuned to give the desired value for following Albrecht and Skordis 
(2000), we have examined a model with A = 0.01,7? = 34.8,7 = 2, and A = 8, and with initial 
conditions chosen so as to fix Q m o = 0.4 today. The evolution of n for this model is also shown 
in Figure 1. It exhibits a sharp transition from n ~ 3.5 at z > 10 to n « 0 at z < 3, a much 
more dramatic change than that in the power-law scalar field models. [Exponential potentials can 
also be made to work in models in which the scalar field is coupled to matter (Amendola 2000); 
however, we confine our attention in this paper to minimally-coupled fields]. 

The dynamical significance of the quintessence component is quantified by the density pa¬ 
rameter Figure 2 shows the evolution of 14^(z) for a cosmological constant (solid curve) 

and the five quintessence models discussed above. The n = 0 (cosmological constant), n = 1, 
and n = 2 cases are quite distinct, as one would expect from their differing p ( p(z). However, the 
V(cj)) = cf)~ l case closely parallels the constant n = 1 case, and the = 0 -6 case likewise tracks 
the model with constant n = 2. The Albrecht-Skordis 14^ is nearly indistinguishable from that of a 
cosmological constant except at high redshift, where the change in n makes a small but noticeable 
difference. 


3. The Observables 


Our starting point is the Friedmann equation in the form of equation (1). As components 
we consider matter with = 0.4, radiation with f4. ri o = 9.8 x 10~ 5 , and quintessence with 

^<70 = 1 — I4m,o — f4 r o ~ 0.6. We compute the ratio p,p(z )/from equation (2) for constant¬ 
'll models [thus obtaining equation (4) for H(z)\ or by computing the evolution of 4> from the 
dynamical equation (5) for the negative power-law or Albrecht-Skordis models. 


The Friedmann equation directly determines the behavior of our first observable, the Hubble 
parameter H(z). We compute other observables given H(z) via the standard treatments in, e.g., 
Peebles (1980, 1993), Kolb & Turner (1990), or Hogg (1999). The age of the universe at redshift z 
is 


t{z) 


dz' 


(1 + z')H{z') 


( 10 ) 


The angular diameter distance dA(z), which is the ratio of the comoving size of an object to its 
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Z 

Fig. 2.— Q'p as a function of redshift z for the indicated quintessence models. The solid curve is 
(z) for a cosmological constant. 


angular size in radians, is 


d A (z) = 


dz' 


H{z>) 


( 11 ) 


(1 + z ) 

Other distances, e.g., those that affect the SN la Hubble diagram or gravitational lensing predic¬ 
tions, are related to d A by powers of (1 + z); the bolometric luminosity distance, for example, is 
d L (z) = (1 + z) 2 d A (z). Since these factors are independent of the cosmological model, a measure¬ 
ment of any of these distances determines all of them to the same fractional accuracy, so we take 
d A (z) as our representative observable for all distance measures. 


The linear growth factor D\ is defined by the relation 

<5 (1) (x,f) = $(x)Z?i(t), (12) 

where S^(pc,t) is the first-order density perturbation. We choose the normalization D\(z = 0) = 1, 
so that D\(z) = S^(z)/S^ (0) gives the linear growth of perturbations between redshift z and 
redshift 0. Then D\ is the growing-mode solution to the differential equation 

D\ + 2H(z)Di — —S2 m) oHg(l + z) 3 D± = 0. (13) 

For fixed Di (z) is a function only of H(z), so it is again determined by the Friedmann 

equation. We solve this equation for D\{z) with a standard Runge-Kutta integration method. In 
the pure cosmological constant case, a closed form expression for D\(z) is 


Di(z) 


H{z) I 

r0 ° dz'{l + z') 

' r°° dz'(i + z')' 

H 0 J, 

* H 3 (z') 

[Jo H 3 (z') \ 


-l 


(14) 
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where the factor in brackets enforces our normalization definition (Eisenstein 1997, based on Heath 
1977). Unfortunately, this expression generalizes only to the case n = 2, but it does illustrate 
that the linear growth factor weights the expansion history in a different way than does the age or 
distance. For constant w, a solution for D\ (z) can be found in terms of hyper geometric functions 
(Silveira & Waga 1994). 

In addition to these four observables, we consider the specific combination h(z ) that is probed 
by the AP geometrical test. AP pointed out that, while tests using cIa(z) or cIl(z) can be affected 
by evolution in the sizes of “standard rulers” or the luminosities of “standard candles,” one can 
measure the ratio of redshift separation distance to angular separation distance assuming only 
that the structures under investigation are isotropic. Recent implementations of this idea consider, 
instead of the idealized spherical clusters discussed by AP, the statistical pattern of clustering 
traced by quasars (Phillips 1994; Ballinger, Peacock, & Heavens 1996; Matsubara & Suto 1996; 
Popowski et al. 1998; Outram et al. 2001; Calvao, de Mello Neto, & Waga 2001), galaxies (Ryden 
1995; Nakamura, Matsubara, & Suto 1998; Nair 1999; Matsubara & Szalay 2001), or the Lya forest 
(Hui, Stebbins, & Buries 1999; McDonald & Miralda-Escude 1999; McDonald 2001). Adopting the 
notation of Phillips (1994) and Popowski et al. (1998), we define 


h 


Az 

TKe 


(15) 


as the ratio of redshift separation to a “redshift arc length” for equal tangential and line-of-sight 
separations in physical coordinates, assuming A z <C z. For a fixed physical separation, Az is 
proportional to H(z) and (A0) _1 is proportional to cIa(z), so h(z) is proportional to their product, 


h(z) = ^-^-H(z)d A (z). (16) 

cz 

We will refer to h(z) as the AP parameter. 

What are the prospects for measuring these observables in the next 5 — 10 years? Our remarks 
here will be qualitative and somewhat speculative, but it is useful to approach the predictions of 
§4 with some sense of what may be achieved by different methods. 


The prospects for distance measurements are the clearest and most well studied. The rms 
scatter of the relation between peak luminosity and light curve shape for SN la is only ~ 10% 
(Phillips 1994; Riess, Press, & Kirshner 1996), so each well observed supernova allows a distance 
estimate with a lcr statistical uncertainty ~ 5%. Current samples (Hamuy et al. 1996; Riess et al. 
1998; Perlmutter et al. 1999) are concentrated at z ~ 0 and z ~ 0.5 — 0.8, but the approach can 
be extended to z ~ 1.7 if the SNAP satellite is built. With samples of hundreds or even thousands 
of supernovae, the statistical errors will become very small indeed, and the measurement accuracy 
is likely to be limited by systematic uncertainties such as dust extinction, possible evolution of 
the progenitor population, and stability of photometric calibration over a wide dynamic range. A 
precision ~ 1% to z ~ 1.7 seems plausibly achievable, and perhaps even unduly pessimistic. 
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There are numerous other ways to measure the distance-redshift relation. Measurements of the 
Sunyaev-Zel’dovich decrement and X-ray properties of clusters can be combined to yield the angular 
diameter distance (Birkinshaw 1999; Molnar, Birkinshaw, & Mushotzky 2002). This method can 
provide an entirely independent check on SN la results, and its limiting precision depends on 
the size of well-observed cluster samples and the accuracy with which internal cluster properties 
(particularly substructure) can be understood. The angular diameter distance can also be measured 
by identifying a characteristic scale (such as the curvature scale of the CDM power spectrum) in 
the angular clustering of distant clusters (Cooray et al. 2001), or even by using the amplitude 
of this angular clustering in comparison to theoretical predictions. Roukema & Marnon (2000) 
have already applied a similar approach to a putative feature at ~ 130 h~ l Mpc in the quasar 
power spectrum. Galaxy counts in a deep redshift survey depend on the volume element dV oc 
d\(z)H~ 1 (z), providing yet another way to constrain distances. The chief uncertainty in this 
approach is evolution of the galaxy population, but Newman &; Davis (2000, 2001) argue that 
this can be controlled in the DEEP redshift survey by measurement of galaxy circular velocities, 
allowing useful constraints on w. 

At redshifts z > 2 all of these methods become difficult, but the Lya forest offers an alternative 
probe out to z ~ 4. The predicted correlation of flux along lines of sight to quasar pairs depends 
on dyi(z), and measurements of this correlation will improve as more close pairs are discovered and 
studied. Such estimates of dyi(z) would be somewhat model dependent, but the statistics of flux 
along individual lines of sight can provide detailed checks of the assumed model. This method 
has not been investigated in any detail (though McDonald 2001 presents relevant results), so it is 
hard to know what precision can be achieved, perhaps a few percent. Strong gravitational lensing 
statistics also test the equation of state through their dependence on distance at various redshifts 
(Cooray & Huterer 1999), and constraints can also be obtained by measuring the source redshifts in 
well understood lens systems (Yamamoto et al. 2001). Finally, the first acoustic peak in the CMB 
power spectrum gives a high-precision measurement of the angular diameter distance to the surface 
of last scattering, at z = z r ~ 1100. The uncertainty in this determination is associated with 
the uncertainty in the parameters that determine the sound horizon at z r , which are themselves 
constrained by the CMB power spectrum and other cosmological observations. It is again hard to 
know just what precision will be obtained on dA{z r ) itself, but a percent or better seems plausible. 

Out to redshift z ~ 1, the main observational probe of the growth factor Di(z), at least 
at present, is the mass function of galaxy clusters. Because clusters are rare objects that form 
from 2 — 3cr excursions of the initial Gaussian fluctuation spectrum, their predicted abundance 
is sensitive to the normalization of that spectrum, and thus to the product a$Di(z), where as is 
the mis fluctuation of matter in spheres of radius 8 h^ 1 Mpc at z = 0 (see, e.g., Bahcall, Fan, 
Sz Cen 1997). However, because the cluster mass function is steep, abundances are also sensitive 
to the accuracy and precision of mass determinations (Frenk et al. 1990), and the limitation on 
measurements of D i is likely to be systematic rather than statistical. The combination of X- 
ray, Sunyaev-Zel’dovich, and weak lensing approaches should reduce these systematic uncertainties 
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below current levels. To guess what level of precision is achievable for D\, we note that current 
“2 ct” uncertainties in the fluctuation amplitude as (for specified £l m , o) are ~ 10% (e.g., Eke, Cole, 
& Frenk 1996), though independent estimates can differ by more than this amount even when 
the input data are similar (see, e.g., Seljak 2001). Balancing the difficulties of working at higher 
redshift against the anticipated large improvements in cluster data, it seems reasonable to hope for 
~ 5% precision in D\ out to z ~ 1, possibly better. Recent discussions of the potential of galaxy 
cluster surveys for constraining w include Haiman, Mohr, & Holder (2001), Newman et al. (2001), 
and Weller, Battye, & Kneissl (2001). 

Cosmic shear is another potential probe of D\(z), measuring the amplitude of surface density 
fluctuations (see Huterer 2001 for a discussion in the context of w constraints). Recent measure¬ 
ments already yield a constraint on as (for fixed ll m ,o) that is competitive with determinations from 
the cluster mass function, with remarkably good agreement of independent estimates (see Maoli et 
al. 2001 and references therein). Measurements of shear for samples of foreground and background 
galaxies with different photometric redshifts should allow as and D\(z) to be disentangled. Weak 
lensing will, at the least, provide an independent check on estimates of D±(z) from cluster masses, 
and the ambitious surveys now underway may eventually yield significantly better precision. At 
z > 2, the most promising route to D\(z) is the flux power spectrum of the Lya forest, which 
is related to the underlying matter power spectrum in a fairly straightforward way (Croft et al. 
1998, 1999, 2001; McDonald et al. 2000; Gnedin & Hamilton 2001; Zaldarriaga, Hui, & Tegmark 
2001). Current uncertainties in the rms fluctuation amplitude are ~ 15%, with roughly equal sta¬ 
tistical and systematic uncertainties. The former will decrease with larger samples such as those 
from the Sloan Digital Sky Survey (York et al. 2000), while the latter will decrease with improved 
determinations of the mean Lya flux decrement, improved numerical simulations to calibrate the 
relation between the flux and matter power spectra, and the use of other statistics to test the as¬ 
sumptions that enter these simulations. At this point, it is not clear where systematic uncertainties 
will limit the precision of mass fluctuation measurements from the Lya forest, but 5 — 10% seems 
a reasonable guess. The other redshift at which we can expect to determine D\(z) is the redshift 
of recombination, from comparing the amplitude of the CMB power spectrum to that of today’s 
matter power spectrum (Doran, Schwindt, & Wetterich 2001). Here measurement precision will be 
high, and the limiting factor is the degeneracy of the fluctuation amplitude with other parameters 
that affect the level of CMB anisotropy. 

We note in passing that the cluster abundance, cosmic shear amplitude, and Lya flux power 
spectrum are not “pure” measurements of Di(z), since the distance-redshift relation affects the first 
two through volume factors and lensing geometry, respectively, and the Hubble parameter affects the 
third because the power spectrum is measured in kms -1 units at the observed redshift. Similarly, 
the angular diameter distance is needed to identify angular scales in the CMB with lengthscales at 
z = 0. However, given the direct dependence of these quantities on the mass fluctuation amplitude, 
it makes sense to describe them primarily as probes of the growth factor. Going from an amplitude 
of fluctuations at redshift z to a value of D\(z) also requires accurate knowledge of the fluctuation 
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amplitude today (i.e., of eg), which we are implicitly assuming will emerge from the tightening web 
of CMB, large scale structure, cluster, and weak lensing constraints. The obtainable precision on D± 
may be higher than the precision in erg itself in the case of a differential evolutionary measurement, 
such as the cluster mass function. 

Lower limits to the age of the universe can be obtained by modeling the stellar populations 
of the oldest galaxies observed at a given redshift. This approach has been used to argue against 
= 1 models, for which the age scales as t(z) = fo(l + z)~ 3 / 2 -, even relative to open models, the 
addition of a cosmological constant makes it substantially easier to understand the red colors and 
high stellar mass-to-light ratios of high-redshift ellipticals (e.g., Peacock et al. 1998; van Dokkurn 
et al. 1998). Lima & Alcaniz (2000) have investigated the usefulness of galaxy ages as a constraint 
on w. Given the uncertainties associated with population synthesis modeling and dust extinction, 
precision of 10% or better in t(z) at high z would seem highly optimistic. However, age constraints 
can provide an upper limit on w that allows a consistency check with other estimates. Exploiting 
this limit requires accurate knowledge of Hq, which sets the overall normalization of timescales. 

The most promising targets for the AP test are quasars (Phillips 1994; Ballinger, Peacock, & 
Heavens 1996; Matsubara & Suto 1996; Popowski et al. 1998; Outram et al. 2001), the Lya forest 
towards quasar pairs (Hui, Stebbins, & Buries 1999; McDonald & Miralda-Escude 1999; McDonald 
2001), and galaxies in the Sloan or 2dF redshift surveys (Ryden 1995; Nakamura, Matsubara, & 
Suto 1998; Matsubara & Szalay 2001). The Lya forest approach is elegant, but McDonald (2001) 
shows that h(z) at z ~ 2 — 4 is insensitive to w, and our results below reinforce this conclusion. 
Instead, h(z) at these redshifts provides a good diagnostic of D m) o (and thus ~ 1 — D mj o), with 
little dependence on w if it is less than —0.5 (McDonald 2001). A precise value of o is needed 
to get useful constraints on w with other tests, as we discuss and illustrate below. Most studies of 
the AP test with quasars or galaxies also focus on o and Da,o rather than w. However, Calvao, 
de Mello Neto, & Waga (2001) have examined constraints on w that could be obtained with the 
2dF quasar redshift survey, with encouraging conclusions. They do not present their results in the 
form of precision on h(z), but their projected sensitivity to w must imply fairly good precision at 
z ~ 0.5 - 1. 

The Hubble parameter H(z) is the observable most directly tied to the Friedmann equation (1). 
One way to measure it is by combining the volume-redshift or AP test with estimates of (1a (z). The 
Lya forest offers a more direct route because the width and separation of features is determined 
largely by Hubble flow (Weinberg et al. 1997). Statistics like the threshold crossing frequency are 
sensitive to the difference between open and flat CDM models because of the difference in H(z) 
(Weinberg et al. 1999a), and measurements of the power spectrum shape can yield characteristic 
scales in kms -1 units at the observed redshift, for comparison with scales measured in hr 1 Mpc at 
z = 0 (Croft et al. 2001). This method of measuring H(z) has not been investigated in any detail, 
so we do not know what precision is attainable; it is likely to be set by the tradeoff between H(z) 
and other parameters that describe the temperature-density relation of the diffuse intergalactic 
medium. It is likely to work better at z ~ 2 — 4 than at lower redshifts, where the observations 
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must be done from space and shock heated gas contributes more to the Ly a forest (Dave et al. 
1999), though even here the separation between features might prove a useful diagnostic of the 
expansion rate. 

At z ^ 1, the skewness of the cosmic shear distribution offers an alternative probe of H(z). Hui 
(1999) discusses the constraints on w that can be obtained by this method, which arise from the 
sensitivity of the predicted skewness to the value of f2 m (z). Since the matter density is necessarily 
Pm(z) = D m; o/9 c ,o(l + z) 3 , the cosmology dependence of £l m (z) comes from the critical density 
p c (z) = 3H 2 (z)/8itG, so in the context of our discussion it makes sense to view weak lensing 
skewness as a measurement of H(z). Jimenez &; Loeb (2001) have proposed yet another route 
to measuring H(z), using the relative ages of galaxy populations at two different redshifts (which 
can be determined more accurately than the absolute ages, since some of the uncertainties in the 
population synthesis models cancel out). The ratio of redshift difference to age difference yields 
dz/dt = —(1 + z)H(z), where the equality uses the definitions (1 + z) = ao/a and H = a/a. Note 
that, while the Lya forest and weak lensing methods effectively measure the ratio H(z)/Hq, the 
age difference method gives H(z) in physical units. Jimenez & Loeb (2001) argue that percent-level 
precision in H(z) is achievable, in which case the uncertainty in the ratio (which is the quantity 
sensitive to the equation of state) is likely to be dominated by the uncertainty in Hq itself. 

There is significant degeneracy between the value of Cl m o and the value of w, since either 
lower fl )rit o or lower w leads to greater acceleration. We assume that improving CMB and large 
scale structure measurements will allow a precise determination of Q m ,o in the next few years, 
independent of measurements of dA, Di, t, H, and h, so that the power of these constraints 
can be brought to bear entirely on the equation of state. We will consider the impact of a 0.05 
uncertainty in the value of D mj o, and it is not obvious whether this assumption is optimistic or 
pessimistic. Apart from the determination of the only role that we ascribe to the CMB is 

the measurement of dA and D\ at z ~ 1100. It may be that CMB data can also yield constraints 
on the expansion history and D\ at lower redshifts, via the integrated Sachs-Wolfe (1967) effect or 
lensing of anisotropies (see, e.g., Seljak 1996), but we do not know just what these constraints will 
be. If the quintessence field is inhomogeneous, it will contribute to large-angle CMB anisotropy 
(Caldwell, Dave, & Steinhardt 1998), allowing a probe of the dark energy independent of the ones 
considered here, which are all based on the expansion history. 

By focusing on this specific set of observables, we do not wish to imply that this is necessarily 
the ground on which theory and observation will be compared. Presumably the constraints on the 
equation of state from, say, weak lensing will be derived in terms of the weak lensing observables 
themselves, without first extracting constraints on D \, dA, and H at various redshifts. However, 
in trying to understand the potential power of combining different observational approaches, it is 
helpful to think in terms of the fundamental quantities that they can measure. In particular, two 
models that predict indistinguishable results for dA(z), D\(z), t(z), and H(z) cannot be discrim¬ 
inated by any combination of observations that depend only on these quantities. Our focus on 
fundamental observables is also a helpful way of estimating the level of precision needed for some 



observational strategy to make a useful contribution to constraining the equation of state and its 
history. 


4. Dependence of the Observables on the Equation of State 

Figure 3 presents our basic results for the quintessence models discussed in §2. Each panel 
shows the evolution of one of the five observables, cU, Di, t, H, or h, out to redshift z = 10 3 ~ z r . 
Open triangles and squares represent constant-w models with n = 1 and n = 2 (w = — 2/3, —1/3), 
respectively. Filled triangles and squares represent V = (j)~ l and V = models, and stars 
represent the Albrecht-Skordis model with the parameters stated in §2. In all cases, the points are 
computed assuming a flat universe, £l mj o = 0.4, H r .o = 9.8 x 10~ 5 , and = 1 — Q mj o — fl r ,o- 
Furthermore, we normalize the value of each observable to the value predicted by a pure-A model 
(n = 0, w = — 1) at the corresponding redshift. Thus, for example, the open triangle at z = 0.5, 
<1a{z)/(Ia{.z )a = 0.96 implies that a precision of 4% (at the desired confidence level) is sufficient 
to distinguish an n = 1 model from a A-model using the angular diameter distance at z = 0.5, if 
£l m o is known perfectly. The ends of the error bars on each point show results for models with the 
same equation of state but Q m fl = 0.35 and 0.45, to illustrate the impact of uncertainty in 
(Ratios are still computed relative to an Q m> q = 0.4 A-model.) If the error bar on an observable 
overlaps a ratio of 1.0, then even a perfect measurement of that observable at that redshift will not 
distinguish the model from a A-model unless fl mj o is known to better than 0.05. 

For a given precision, the sensitivity of different observables peaks at different redshifts. The 
Hubble parameter sensitivity peaks at z ~ 1 — 2, when the ratio of fl^(z) values in different mod¬ 
els is large and the quintessence energy density is still large enough to be dynamically important 
(see Figure 2; roughly speaking, it is the absolute difference in in this Figure that matters for 
differences in H). The sensitivity of cIa(z) remains fairly flat, since even at high z the distance 
“remembers” the behavior of H(z) at low redshifts (see eq. 11). The age of the universe, by con- 
strast, depends only on the Hubble parameter at redshifts higher than z (eq. 10), so the sensitivity 
of t(z) continues to increase almost all the way down to z = 0. Note that the sensitivity of an 
observable to the value of n, displayed in Figure 3, may be quite different from the sensitivity of 
that observable’s z—derivative, which often peaks at lower redshift. 

The behavior of h(z) is governed by the competing effects of H(z) and cIa(z) — from equa¬ 
tions (16) and (11), one can see that h(z ) is proportional to the product of H(z) and the average 
value of H^ 1 (z) at lower redshifts. At z ~ 3, h(z) is very insensitive to w, as pointed out by Mc¬ 
Donald (2001), who emphasizes that this independence of w makes the AP test at this redshift an 
especially good diagnostic of D m> o- The sensitivity to w peaks at z ~ 0.5, making quasar clustering 
better than the Lya forest as a probe of the equation of state per se. The Sloan survey’s luminous 
red galaxy sample (Eisenstein et al. 2001) might also be useful for this application (Matsubara & 
Szalay 2001). The sensitivity of h(z) to w grows again at z > 5, but the prospects for applying the 
AP test at these redshifts seem very slim. 
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Fig. 3.— Di(z), t(z), (1a(z), H(z), and h(z) versus z for the constant n = 1,2 models, the fixed 
potential models V = and the Albrecht-Skordis potential. The central value represents 

Q m p = 0.4 and the error bars are at = 0.35, 0.45. Quantities are normalized to the A model 
with o = 0.4. A small horizontal offset has been added to the points to allow them to be 
distinguished. The fl m o = 0.45 end of the error bar is usually the end further from a ratio of unity 
(or, for the A-S model, the end closer to the n = 1 points), except for D\(z) and the high redshift 
(. h < 1) regime of h(z), where the il m fl = 0.45 end is closer to unity. 
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The sensitivity of D\ (z) increases with increasing z, since the growth factor depends only on 
clustering between redshift 2 and redshift zero, and models with lower n have larger Q m (z) = 
1 — Q^z) — Q r (z) at all z > 0. For n = 1, the sensitivity levels out at z > 3 as quintessence 
becomes dynamically unimportant, but for n = 2 the value of Cl^z) is non-negligible even at fairly 
high redshift (Figure 2). 

Figure 3 shows that ~ 10% measurements of any of these observables near their redshift of 
peak sensitivity can discriminate an n = 2 model from a A-model, or from an n = 1 model. Since 
the expansion of an n = 2 model with fi m ,0 > 0 is always decelerating, it is not surprising that this 
model is fairly easy to distinguish from a A-model with significant acceleration at low redshift. SN 
la measurements already rule out this value of n (Garnavich et al. 1998), and Figure 3 implies that 
other observations are within reach of confirming this result independently. Distinguishing an n = 1 
model from a A-model is much harder, typically requiring measurement precision of a few percent 
or better, and independent precision on 0 that is not much worse than the 0.05 represented by 
our error bars. Nonetheless, this level of discrimination is clearly within reach of the improving SN 
la measurements, and the discussion in §3 suggests that several other methods have a realistic hope 
of reaching the necessary precision on a timescale of several years. The sensitivity of our observables 
to w would be slightly greater if we adopted fl mi o = 0.3 as our central value. Furthermore, the 
precision that might be obtained from measurement of these observables at multiple redshifts can 
be higher if the errors in the separate measurements are uncorrelated; however, if the source of 
uncertainty is systematic, it may produce correlated errors at different redshifts. 

Unfortunately, distinguishing any of the time-varying w models that we have considered from 
the closest constant-^ model looks all but impossible. The V = q U 1 model tracks the n = 1 model 
almost perfectly, and the V = model tracks the n = 2 model with similar faithfulness. The close 
match of these models is unsurprising given the plots of n(z) and in Figures 1 and 2. The 

predicted differences between the power-law potential models and the constant-n models, and the 
redshift-dependence of these differences, have the sign one would expect from Figure 1; the problem 
is simply that the time-dependence of the equation of state predicted by these models is extremely 
weak. Similarly, the Albrecht-Skordis model is virtually indistinguishable from a pure-A model 
because it has n ~ 0 at all redshifts where quintessence is dynamically important, even though it 
has a very different n at z > 3. The one potential distinguishing feature of the Albrecht-Skordis 
model is the value of D\ at recombination, which is about 7.5% larger than that of a A-model. 
This level of precision is plausibly within reach of future observations. The distinguishability of 
the Albrecht-Skordis model would increase if the equation-of-state transition were shifted towards 
lower redshift, and vice versa. 

Since we have so little empirical information about the nature of dark energy, there is no reason 
to think that models presently in the literature exhaust the possibilities for the time-dependence 
of the equation of state. We have therefore constructed a set of “toy” models that exhibit a wider 
range of behavior, so that we can better understand the ability of observations to detect time- 
variation if it is present. For these models, we assume that the redshift-dependence of is a 
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broken power-law of the form 


P<t>( z ) 


(1 + z c ) n ~ m (l + z) n 
{l + Z C + z) n - m 


0 ), 


(17) 


where z c is the critical redshift near which the scaling behavior of p$ changes from (1 + z) m for 
z z c to (1 + z) n for z <C z c . The model therefore switches from an early time const ant-m case to 
a late time constant-n case. This class of models allows us to examine the effects of a more extreme 
change in the power-law scaling of p^ than is exhibited by the tracker models examined earlier. 
(For the cases m = n or z c = 0, this model reduces exactly to the constant-n models discussed 
earlier, while the m = 3, n = 0 case resembles the Albrecht-Skordis model). Note that a similar, 
but somewhat different toy model was examined by Huterer and Turner (2001), who looked at 
models in which w was taken to be constant in discrete redshift bins. 


Figure 4 shows fl<f>(z) for a variety of broken power-law models. The late time behavior (n) 
has been fixed at n = 1, and early time behavior (m) has been taken to be m = 0 or m = 2, 
for two different values of the critical redshift (z c = 1 and z c = 3). For comparison, we have also 
included three constant-n models (m = n), namely m = n = 0,1,2. As expected, fl^( z ) in our 
broken power-law models deviates from its behavior in the constant-n models to a much greater 
extent than is the case for the power-law potentials in Figure 2. 

In Figures 5—7 we examine our five observables for the broken power-law cases ?n, n = 0,1,2 
and z c = 1,3. Each page of graphs shows models with a different late time behavior (different 
value of n), and the observables on each page are normalized to the corresponding constant-n case: 
n = 0 in Figure 5, n = 1 in Figure 6, and n = 2 in Figure 7. The deviation from a ratio of unity 
in each case shows the observational effect of the break in scaling behavior. This deviation is quite 
significant in many cases, often more than 10% for H , d,A or t, and up to almost 30% for D\. The 
broken power-law model is therefore clearly distinguishable from the constant-?! model that has the 
same value of n at z = 0. 


However, if we choose a constant-n model that is matched to the “average” behavior of the 
broken power-law model instead of the z = 0 value of n, then this distinguishability vanishes. The 
curves in Figures 5—7 show, for each broken power-law case, the predictions of a constant-n model 
selected to produce the same value of the Hubble parameter H(z) at z = 1. This matched constant- 
n model predicts nearly the same values for all observables at every redshift as the corresponding 
broken power-law model. By design, our toy model has a large change in n (and thus w) at low 
redshift, but this time-variation cannot be detected unless the observables can be measured to 
extremely high precision. Making the transition redshift z c higher or lower does not make the 
time-dependence easier to detect; it just changes the effective average value of n, so the constant-n 
model that matches the broken power-law model is different from before. For example, we have 
constructed models with z c = 0.5, and our results are quite similar, except that the corresponding 
constant-n model has a value of n which is closer to the low-redshift exponent in the broken power- 
law model. A sharper transition would be easier to detect, but our models already change the 
energy scaling index by order unity over a redshift interval Az ~ z c [and thus a time interval 
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Fig. 4.— as a function of redshift z for the indicated broken power-law models, defined by 
eq. (17). The curves show three constant-?} models for reference. 

At ~ H~ 1 (z)], and a much faster transition seems physically unlikely, although such models can 
be constructed (see, e.g, Weller and Albrecht 2001). 

While our results suggest that generic variations of the equation of state with redshift are 
essentially undetectable, there is one exception to this rule that is worthy of note. Our models with 
high-redshift index rri = 2 predict values of D\{z) at z = 10 3 that differ by a few percent from those 
of the constant-?? model that matches the low-redshift observables. This level of precision might 
plausibly be achieved with comparisons of CMB anisotropy to local clustering, though substantial 
improvements in observational data would be required. [The values of t(z), H(z), and h(z) also 
show percent-level deviations from the matched constant-?} model at z > 3, but we see no plausible 
routes to attaining the necessary precision for these quantities.] The behavior of the m = 2 models 
is reminiscent of the Albrecht-Skordis model, and the cause is similar: the energy density of the 
quintessence component remains a non-negligible fraction of the critical density out to fairly high 
redshift, so the gravitational growth of matter clustering is correspondingly slower. 



5. Conclusions 

Our results both confirm and extend earlier work on this subject. For any of the five observable 
quantities considered here — the angular-diameter distance cIa(z), the Hubble parameter H(z), 
the age of the universe t(z ), the linear growth factor D\(z), or the Alcock-Paczynski parameter 
H(z)cIa(z) — measurement with ~10% precision near the observable’s redshift of peak sensitivity 
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Fig. 5.— D\(z), t(z), cIa(z), H(z), and h(z) versus z for four broken power-law models, all having 
the same late-time behavior, n = 0. The points are normalized to the value of the given observable 
in the n = 0 model at the same redshift. For each set of points, a constant-n model that has the 
same value of H(z) at z = 1 is shown. 
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would be sufficient to distinguish an n = 2 (w = —1/3) model from a pure cosmological constant, 
even if fl m) o were known only to an accuracy of ±0.05. Although this value of w is already ruled 
out by the SN la measurements, our results suggest that other observations may soon be able 
to independently confirm the result. Distinguishing an n = 1 (w = —2/3) model from a pure 
cosmological constant is much harder, demanding measurement precision of a few percent near the 
redshift of peak sensitivity, along with a determination of D mj o to within ±0.05. Although this 
level of precision is currently unavailable, it seems clearly within reach of improving SN la data, 
and it is likely to be achieved by one or more of the other observational methods discussed in §3. 
Thus, while SN la surveys may provide the first precise determination of w, a collection of other 
observations seems likely to provide confirmation (or refutation!) of the measurement within a few 
years. 

The sensitivity of the observables to the value of n depends on redshift in different ways, 
reflecting the links between these quantities and the expansion history. The age t(z ) depends 
only on expansion at redshifts greater than z, so its sensitivity to n decreases monotonically with 
increasing z. The linear growth factor, on the other hand, depends on clustering from redshift z 
to redshift zero, so the sensitivity of D\(z) increases monotonically with z. The Hubble parameter 
H(z ) is most sensitive at z ~ 1 — 2, when is substantially different from its present-day value but 
not so small that quintessence is dynamically unimportant. The sensitivity of the angular diameter 
distance is fairly flat over a wide range of redshift. The sensitivity of the AP parameter is governed 
by competing effects of H(z) and cIa(z), which cancel each other at z ~ 3. 

Because of their different connections to the expansion history, we hoped at the outset of this 
investigation that these observables would provide complementary information about the history 
of the equation of state, allowing a combination of measurements to detect a time-variation of w 
that could not be found by any one method on its own. Unfortunately, we find that the level of 
complementarity is too weak to be useful in practice: models that make indistinguishable predic¬ 
tions for one observable generally make indistinguishable predictions for all of them. Of course, it is 
valuable to confirm an important result like a measurement of w by independent methods, to check 
for systematic errors or a breakdown of the assumptions implicit in each approach. Also, different 
observables can provide complementary information about Cl m o, precise knowledge of which is es¬ 
sential if one hopes to constrain w. However, once Q m o is known, the constraints on the equation of 
state and its history will be dominated by the single highest precision measurement; adding lower 
precision measurements of other observables will give little additional purchase. 

We find, furthermore, that none of the observables holds much promise for distinguishing a 
quintessence model with a time-dependent equation of state from an appropriately chosen constant- 
n model, even if one is highly optimistic about the achievable precision and assumes perfect, in¬ 
dependent knowledge of Tracker models with V{(j)) oc and V{4>) oc </> -6 are effectively 

identical to models with constant n = 1 and n = 2, respectively. Models with an Albrecht-Skordis 
potential cannot be distinguished from a pure-A model, except, perhaps, by a measurement of the 
growth factor at recombination from CMB anisotropy (a point we will return to shortly). The fun- 
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damental difficulty is that, in any observationally viable model, quintessence becomes dynamically 
important only at low redshift, so it affords little purchase for measuring redshift dependence of 
its equation of state. Furthermore, as Figures 5—7 demonstrate, even a substantial transition in 
n at low redshifts is very difficult to detect, since the value of n at 2 ; = 0 is not known a priori, 
and time-variation must therefore be judged relative to the constant-n model that best mimics the 
time-variable model. Our broken power-law models have substantial low redshift transitions by 
design, but there is usually a constant-n model that predicts the same values of all observables 
to better than 1% at all observationally accessible redshifts. We conclude that detecting time- 
variation of the equation of state requires truly extraordinary precision, unless this variation occurs 
on a timescale much shorter than the Hubble time, which is possible but seems physically unlikely. 
Sub-percent precision may be achievable by some methods (SN la observations look to us like the 
best hope), but it requires controlling systematic uncertainties, especially those that are correlated 
among different redshift bins, very tightly. 

Our conclusions in this regard agree with those of Maor, Brustein, & Steinhardt (2001), who 
found that accurate measurements of the luminosity distance alone would be insufficient to de¬ 
termine the form of w(z ) for the dark matter energy component. Wang & Garnavich (2001) and 
Tegmark (2001) showed that SN la measurements should be able to detect time-variation in the 
energy density p,p(z), but this only means demonstrating that n > 0 (w > —1); we agree that 
a significant departure from n = 0 should be detectable, but detecting time-variation of n is far 
more challenging. In a similar vein, despite fairly optimistic assumptions about the prospects for 
the SNAP satellite, Huterer & Turner (2001) find that error bars on the time-derivative of w are 
quite large, and degrade considerably with uncertainty in Yamamoto et al. (2001) suggest 
that the form of the dark energy equation of state might be determined by studying strong grav¬ 
itational lensing systems, but their results indicate that detecting time-variation is possible only 
with extremely high precision measurements of the lensing systems, and then only if fi mj o is known 
precisely. The principal significance of our results, relative to these earlier papers, is that they ap¬ 
ply to all proposed observable tests based on the cosmic expansion history, since these tests always 
measure some combination of H(z), cIa{z), t(z), or D\(z). 

Our investigation shows that there is one generic form of time-variation in the equation of state 
that might be observationally detectable. Constant-n models with n > 2 (u> > — 1/3) are ruled out 
by current data, but a time-variable model could have n > 2 at high redshift and a transition to 
low n at low redshift when quintessence becomes the dominant energy component. The Albrecht- 
Skordis model displays just this behavior, since the quintessence roughly tracks the matter energy 
density (n « 3) along the exponential part of V((j>) but changes its equation of state (to w « — 1, 
n ~ 0) when it reaches the potential minimum. If n > 2 down to some fairly low redshift, then 
the dynamical effects of quintessence are non-negligible (though small) over a fair fraction of the 
post-recombination expansion history, and they slow the progress of matter clustering. The result is 
a slight (few percent) mismatch between the value of D\ at 2 = z r ~ 1100 and the value expected 
for a constant-n model that matches the low redshift data; in observational terms, the level of 
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CMB anisotropy would be a few percent higher than anticipated. Doran, Schwindt, & Wetterich 
(2001) emphasize a similar point and discuss the relation between CMB anisotropy and erg in detail. 
Detecting even this type of time-variation will be very challenging, requiring a precise determination 
of the effective low-redshift value of n, precise determinations of the present-day amplitude of matter 
clustering and Q m o, and the demonstration that any excess CMB anisotropy does not arise from 
other sources, such as tensor fluctuations, secondary anisotropies, or contaminating foregrounds. 

The discovery of dark energy is an extraordinary cosmological achievement, one that could 
happen only in the era of “precision cosmology.” If the equation of state of this dark energy is 
substantially different from p = —p, or if it has been different in the recent past, then that departure 
should be detected independently by several of the ambitious observational efforts currently planned 
or underway. A precise (~ ±0.1) measurement of the low-redshift value of w would be another 
extraordinary achievement, ruling out many models for the origin of dark energy and tightening the 
parameter space of others. However, the information provided by different observable probes of the 
cosmic expansion history, or by the same probe at different redshifts, is mostly redundant rather 
than complementary, once D mj o has been determined to high precision. As a result, the next step 
of detecting time-variation in the cosmic equation of state is likely to prove extremely difficult. If 
we are lucky, then the dark energy has the kind of dynamical significance at high redshift or sudden 
transition at low redshift that produces an observationally accessible signature, though reading that 
signature will still require a combination of several cosmological measurements of unprecedented 
precision. If we are not so fortunate, then the observable effects of the dark energy will, for the 
foreseeable future, provide only two numbers with which to describe it, the current energy density 
and an effective low-redshift value of w (or some equivalent pair of parameters). Until a physical 
model comes along that accounts for these two numbers in a natural way without adjustable inputs, 
the true nature of the dark energy component is likely to remain mysterious. 
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